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Abstract 

We analyze current-voltage characteristics taken on Au-doped indium-oxide films. By fitting a scaling function to the data, 
we extract the electron-phonon scattering rate as function of temperature, which yields a quadratic dependence of the electron- 
phonon scattering rate on temperature from IK down to 0.28K. The origin of this enhanced electron-phonon scattering rate is 
ascribed to the mechanism proposed by Sergeev and Mitin. 

PACS numbers: 72.15.Lh,72.10.-d, 73.61.-r, 63.20.Kr, 72.10.Di 



(N 



S3 

S3 
i 

C/3 



I 

S3 
O 

o 



> 
in 

OS 

in 
in 
o 
in 
o 

c3 



i 

S3 
O 

o 



X 



I. INTRODUCTION 

Energy relaxation processes play an important role 
in the low temperature transport properties of diffusive 
metals, alloys and semiconductors. In particular, they 
determine the maximum electric field F allowed if the 
system is to be measured under near-equilibrium condi- 
tions at a given temperature. The condition for that 
may be expressed as eFL e <C fc.eT where the energy 
relaxation length L e is the length over which the elec- 
tron diffuses under F before the energy gained from it 
is dissipated into the thermal bath. In particular, this 
issue is relevant for all aspects of quantum transport 
such as corrections to the conductivity due to interfer- 
ence and electron-electron interactions. For a system of 
size L ^> L e , energy relaxation processes are usually con- 
trolled by electron-phonon inelastic scattering, and in the 
following we shall assume that L e is dominated by the 
electron-phonon diffusion length L ep . 

In clean samples the electron-phonon scattering mech- 
anism is well understood, and the scattering rate, r~ p l , 
is known to be proportional to T 3 (where T denotes the 
temperature). In dirty systems however, where the elas- 
tic mean free path of the electron is smaller than the 
phonon thermal length, the situation is more compli- 
cated. Schmidi showed that in this case T~ p is sup- 
pressed and becomes proportional to T 4 , in accordance 
with Pipard's ineffectiveness condition 2 .. His model as- 
sumed that the impurities are anchored to the lattice, 
and the scattering rate was calculated by moving into 
a reference frame which follows the lattice vibrations. 
Riezer and Sergeev obtained the same result using the 
laboratory frame of reference 3 -. A scattering rate propor- 
tional to T 4 has indeed been observed in a number of 
experiments^^. However, a T 3 law was frequently ob- 
served even in systems that were presumably well into 
the dirty limit regimeZ&Siii. Moreover, quite a few ob- 
servations of a T 2 scattering law were reported in other 
experimentsii*i2iiiii4. The latter experimental results 
triggered further studies with the aim of understand- 
ing better the electron-phonon scattering mechanism in 
disordered metals. In particular, to obtain an electron- 
phonon scattering rate proportional to T 2 (rather than 
the "ineffective" T 4 law), Belitz and Wybournei^ as- 



sumed a strong phonon damping, while Jan Wu and 
Weii& included effects associated with the discrete lattice 
structure. Sergeev and Mitin 17 obtained the T 2 behavior 
from a model were impurities are assumed to be fixed, 
namely, impurities which do not follow the lattice vibra- 
tions. They argued that heavy impurities or boundaries 
which move differently from the host lattice produce the 
same effect. 

In this paper we analyze the non-ohmic characteris- 
tic of thin films of In2 03_ K :Au (crystalline indium-oxide 
doped with 2% gold), that were characterized and mea- 
sured as described elsewhere^. We show that in this 
system r^ 1 oc T 2 , and interpret this behavior as a mani- 
festation of the Sergeev Mitin mechanism where the Au 
atoms play the role of the "immobile impurities" . The 
sample resistance is used as a thermometer of the elec- 
tron temperature. The latter is determined by the energy 
balance between the Joule heating, due to the presence of 
electric field, and the heat transfer to the lattice phonons 
governed by T~ p x . 



The data we shall analyze are dynamic resistance R 
versus voltage traces for a typical film at T < IK shown 
in Fig. 1 (left panel). As expected, the deviations from 
ohm's law become more pronounced as T decreases. In 
the right panel of Fig. 1 it is shown that these curves 
taken at different temperatures can be collapsed onto a 
common function: 



AR _ R{F,T) - R(0,T) 
~R = R{0,T) 
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with p = 2. The scaling of the electric field as a power 
of the temperature, F ~ Ti" 1 " 1 , and its relation to the 
electron-phonon time T~ p (T) oc T p , has been already rec- 
ognized by Anderson, Abrahams and RamakrishnaniS, 
and later by Arai 2 ^. In the next section we shall cal- 
culate the function ATZ p , and clarify its relation to the 
electron-phonon relaxation time. We shall consider the 
dependence on the sample length and compare with fur- 
ther experimental results in section 3. 
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FIG. 1: The non-ohmic characteristic of 200Athick ft^Oa-^Au sample with length L=3500/im and width W=lmm. Left 
panel: The differential resistance (defined in Eq. Q as function of the voltage. Right panel: The same data plotted as function 
of F(= V/L) normalized by T 2 . Note the near perfect data collapse as well as the fit to formula 12011 represented by the 
continuous line. 



II. THEORY 



The essence of the picture developed below is the as- 
signment of an effective temperature for a given field 
and bath-temperature. The electrons are accelerated 
by the electric field, and the collisions with other elec- 
trons and phonons result in a local equilibrium distri- 
bution characterized by an effective temperature T e ff. 
The latter differs from the bath temperature T, and de- 
pends on the electric field. We shall make the association 
R(F, T) ~ i?(0, T eff ). Hence knowing T eff {F, T) and the 
form of the near-equilibrium R(0, T) yields the desired 
R(F,T) from which we deduce the scaling function (Q. 

We begin by considering the Boltzmann equation for 
the electron distribution function /: 



d£ 
dt 



dr 



eF 



dp 



= I[fl 



(2) 



Here r is the position, v is the velocity, p denotes the mo- 
mentum, F the electric field, and I[f] = Ii m [f] + Iee[f] + 
I ep [/, N] represents the collision integrals due to impurity 
scattering, electron-electron, and electron-phonon inter- 
actions. The latter depends also on the phonon distribu- 
tion function, assumed to be the equilibrium distribution 
function denoted by N. 

Following Nagae^*2i and Kozub and Rudin 22 , we look 
for a steady state solution of the form: 



/ = f(r,h,e- eF ■ r) 



(3) 



where h denotes a unit vector in the direction of the mo- 
mentum, and e = e(p) is the energy assumed to depend 



on the absolute value of the momentum p — \p\. Substi- 
tuting |gj in © leads to: 



dr 



ij 



(4) 



Next, we define the symmetric and anti-symmetric parts 
of the distribution function with respect to the momen- 
tum direction = (f(n) ± /(— h)) /2, and construct 
two new equations from (JlJ associated with the addi- 
tion and subtraction of Boltzmann equations for f(n) 
and /(— n). Then assuming that momentum relaxation 
is dominated by scattering from impurities, and that 
electron-electron and electron-phonon interactions are es- 
sentially independent of the momentum direction (i.e., 
I ee [f] ~ /«[/+], and I ep [f,N] ~ I ep [f+,N]), we obtain: 



r df~ 
dr 



ij 



UiTlj) 



dn 3 



where /[/+] ~ I ee [f+] + I ep [f+,N], and 



dr 



(5) 



(6) 



In the simplest approximation, the impurity collision 
term takes the form Ii rn [f~] = —f~/r where r is the 
elastic mean free time. Thus using (JHJ), the antisymmet- 
ric part of the distribution function can be expressed in 
terms of the symmetric part and substituted into JSJ. 
The resulting equation is now averaged over the momen- 
tum directions to give 



-DV 2 f+ =I ee [f+]+I ep {f + ,Nl 



(7) 



2 



where D = TUp/3 is the diffusion constant (vp is the 
Fermi velocity). Here and henceforth we neglect the en- 
ergy dependence of the diffusion constant. 

We wish to solve Eq. for homogeneous quasi-two 
dimensional samples of rectangular shape with voltage 
contacts located at x — ±L/2. Thus / + is independent of 
the transverse coordinate, and the boundary conditions 
assuming ideal contacts are: 



=±L/2 



np 



(8) 



where np(x) — (1 + exp(x))^ 1 is the Fermi distribution 
function, e is the electron energy, ep is the Fermi energy, 
V = FL is the voltage drop across the sample, k B is 
Boltzmann constant, and T is the bath temperature. 

Equation Q is a nonlinear equation for the electron 
distribution function. To make progress, we shall assume 
that the electron-electron diffusion length L ee is much 
smaller than the system size, L, and the energy relax- 
ation length L ep . This should secure an effective local 
thermalization due to the large number of collisions an 
electron experiences. Looking then for a solution which 
describes local equilibrium of the electrons: 



np 



e — ep — eFx 



(9) 



where Ti oc (x) is a local temperature of the electrons. At 
the contacts the local temperature by assumption coin- 
cides with the bath temperature: 



Tic 



4 1 = - 



(10) 



so that the solution @ with the boundary conditions 
(|f 0|l satisfies the requirement (|HJ). Under the assumption 
of local equilibrium the electron-electron collision term 
vanishes and equation J2J) reduces to 



2 f+ 



-D 



d 2 f 
dx 2 



= Iep[f + ,N] 



(11) 



Finally, to extract the local temperature behavior we 
multiply this equation by e and integrate over the en- 
ergy. The resulting equation is: 



D 



(> 



dx 2 



{eFf 



deeI ep {f+,N}. (12) 



To further simplify l|12(l , we take the electron-phonon col- 
lision integral to have the form: 



(i-n.+^rwiv^) 



(LjK(uj 

(WW(^)(a^) + 

(l-/+(e-a;))/+( e ) (n ( 



k B T 



+ (l-/ + (e))/ + (e-^(-^) 



(13) 



where K(u>) = aw p_1 , is a kernel depending on the na- 
ture of the collision between the electrons and phonons, 
while N(x) = (exp(a;) — 1) _1 is the equilibrium phonon 
distribution function. We substitute this expression into 
(|12fl with the approximate local equilibrium form of the 
electron distribution function JJ]J, and integrate over e 
and lu. The resulting equation is an equation for the 
local temperature 



v 2 k 2 B dT 2 oc (x) 
6 9a; 2 

where 



+ (eF) 2 =^+ 2 fc 2 (*)-T»+ 2 



(14) 



r;=(l-2-^ +1 >)(p+l)!C(p + 2)-. (15) 

In understanding the form of the solution of equa- 
tion (|14|) . it is instructive to identify first the relevant 
length scale for this equation. To this end one may lin- 
earize the equation by substituting T 2 oc (x) = T 2 + 8T 2 (x) 
and expanding to linear order in 8T 2 {x). The solution of 
the resulting equation with the boundary conditions I|1U|) 
is 



T 2 oc (x) =T 2 + 



where 



6(eFL ep ) 
k 2 B n 2 



cosh 



cosh 



L 

2L,, 



(16) 



Tt(k B T) I 



(17) 



is essentially the electron-phonon length at equilibrium. 
From (|16[) one can see that L ep sets the distance from 
the contacts over which the temperature profile reaches 
a saturated value. Thus long sample satisfies L 3> L ep , 
and these may be considered to have an essentially space 
independent local temperature which we shall refer to as 
the effective temperature. 

The linearized solution for long samples l|16|) is strictly 
justified when the electric field is weak, i.e. eFL ep < 
ksT. We shall be interested in the limit of strong fields 
where L ep will presumably be smaller than its equilib- 
rium value, the r.h.s of l|17(l . Provided L ^ L ep (T,F) 
it makes sense to assume an electron temperature, which 
is essentially constant throughout the sample. We then 
neglect the space dependent term in equation (|14fl . and 
the effective temperature of the electrons, at any field 
strength, is readily deduced to be: 



T, 



eff 



(eFf 

vk p B +2 



(18) 



As mentioned earlier, once T e ff is known, to find 
the scaling function A7^. p requires only the temperature 
dependence of the resistance, R(T), since R(F,T) ~ 
i?(0, T e ff). At low temperatures, the temperature depen- 
dent terms of the resistance are the weak localization 23 



3 



and the Altshuler-Aronov corrections^. For thin films, 
both of these have a logarithmic behavior, thus 



fl(G,T)~i? (l- 7 ln(T)), 



(19) 



where 7 is a small dimensionless constant, depending on 
the nature of the electron-electron interactions, the spin- 
orbit coupling, and the ratio of the quantum unit resis- 
tance to the Drude resistance of the sample, Rd- Sub- 
stituting R(F,T) ~ R(0,T eff ) and in the definition 
of the scaling function Q , and expanding to the leading 
order in 7 we obtain: 



ATZ. n 



F 



7 



■In 



P 



(eFf 



v(k B T) p+2 



(20) 



Note that this has the scaling form as in Eq. (JTJ above. 

To anticipate the discussion in the next section, it is 
instructive to consider the case of short samples L ep 3> L. 
Here one may neglect the contribution of the electron- 
phonon collision term in Eq. i|12|) and readily obtain the 
solution for Ti oc (x), whose space dependence, in general, 
cannot be ignored: 



T? oc {x) =T 2 + 



3(eF) 5 



4 



(21) 



At the temperature range where the dephasing length 
is much smaller than the system size one may view the 
sample as a set of classical resistors connected in series. 
Therefore the total resistance can be approximated by 
the sum: R{F,T) ~ J2j R j> where R j = R(Tioc(xj)) is 
the resistance of the j-th segment (of size of the dephas- 
ing length), centered at the point Xj. Thus, assuming 
homogeneous sample, the experimentally measured re- 
sistivity is essentially an average over the position. From 
this average one immediately obtains the scaling function 
of short samples: 



where 



X 



\ 1 .anii 1 ( i 1 - I 



4tt 2 (k B T 



(22) 



(23) 



3 \ eV 

and V — FL is the voltage drop along the sample. 

III. DATA ANALYSIS AND DISCUSSION 

Comparing the resistance curves shown in Fig. 1 with 
Eq. I|20l) . one notes that a scaling form that leads to 
the data collapse occurs for p = 2. This means that 
the energy relaxation time is quadratic with tempera- 
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FIG. 2: Resistance versus temperature for the sample of Fig. 1 
(open squares) and for a h^Oa-^Au sample with length 
L= 80 fim and width W=500/im (full circles). The dashed 
lines are the logarithmic slopes used in defining the value of 
7 (see Eq. JT3J). 



The procedure of extracting the detailed form of T~ p 
or L ep from the experimental results is as follows. First, 
the value of r\ is obtained by fitting the R(F , T) data of 
Fig. 1 to Eq. (|2*U|) . as shown in the right panel of this 
figure. Note that the latter needs the parameter 7 as in- 
put. This is defined by equation (|19l) above and is thus 
obtained from the near-equilibrium R(T) measurement 
performed on the same sample. Such R(T) data and their 
associated 7 are shown in Fig. 2 for two samples. These 
are made from the same batch of a Au-doped IniO-j,- x 
film, they only differ in their lateral dimensions. The 
first is the 3500/_tm sample of Fig. 1. The second sample 
is 80/iTO long. The sheet resistances of the two are within 
1% of each other, yet their logarithmic slopes are some- 
what different. Also, both samples show a systematic de- 
viation from the theoretical ln(T) dependence, a feature 
that seem to occur in some other 2D systems^. Since we 
are mainly interested in the restricted temperature range 
0.28 — IK, this feature may be ignored, and 7 is defined 
by fitting R(T) to a simple ln(T) over the relevant range 
as shown in the figure. The fits yield 7 ~ 0.0098 and 
7 ~ 0.0081 for the long and short samples respectively, 
and these values are used in the subsequent analysis be- 
low. 

An excellent fit to the data in the right panel of 
Fig. 1 can be obtained using Eq. I|20|) with 77 = 3.7 • 
10 55 , h i , 7 = 0.0098 and p = 2. The use of this for- 

Joule* 1 ' 1 

mula, appropriate for the L 3> L ep limit, is justified for 
this sample as can be seen by estimating L ep . Inserting 
the above value of 77 in equation l|17f) . gives L ep ~ 20 fim 



ture: r ep oc j^- oc T , and therefore the electron-phonon at 1K and only L ^ „ 6Qflm at T ^ Q 28K Thus L 
length l|17fl is inversely proportional to the temperature. much smaller than L down to the lowest temperature we 
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FIG. 3: Dependence of the differential resistance on volt- 
age at several temperatures. In203_ x :Au sample with length 
L—80fim and width W=500/im. The noisier data here (as 
compared with those of Fig. 1), is mainly due to the much 
smaller sample size. 



are dealing with here. 

For the 80fj,m sample, on the other hand, a crossover 
to the short-sample regime is realized in the temperature 
interval covered in our experiments. The crossover can 
be seen by studying the data depicted in Figs. 3, and 4. 
Fig. 3 shows the raw R(V) data measured at different 
temperatures. Fig. 4 show these data plotted according 
to the short-sample formula (Eq. (|22l) ') in the left panel, 
and according to the long-sample scheme (Eq. I|2(J[) with 
p = 2) in the right panel. The crossover temperature is 
the temperature below which the R(F, T) data can be 
scaled as a function of V/T, and above which it scales as 
F/T 2 . Comparing between the left and the right panels 
of Fig. 4 it is evident that this temperature is approxi- 
mately 0.5K . Thus, with this sample, the consistency of 
our approach can be tested in the two limits. As shown 
in the left panel of Fig. 4, a good fit to the AR/R(0) 
data is obtained using Eq. (|22|l . which involves just the 
parameter 7. Note that the best fit 7 is quite close to the 
7 that one gets from the logarithmic fit to the R(T) data 
of this sample (Fig. 2). The other limit, which conforms 
to Eq. l|2U)l . also yields a reasonable agreement, with the 
same 7 as used above, and with ri = 2.2 • 10 55 -? — ^ — r 
(see right panel of Fig. 4). The quality of the fit here 
is less good than in the 3500/Ltm sample, perhaps due to 
the fact that even at the highest temperature used the 
sample is not really in the long-sample limit. 

Finally, using the value of 77 for this sample in Eq. (|T7Jl . 
one gets L ep ~ 34/im at IK and~ IOO/xto at T !~ 0.28K . 
Since it is plausible to expect that the crossover from the 
short-sample to long-sample regime should occur when 
L ep becomes comparable with L/2, these numbers are 



consistent with our picture. 

We turn now to discuss the physics that underlies 
the r~ p oc T 2 law for the electron-phonon scattering 
rate suggested by our analysis. It is important to note 
that the enhanced electron-phonon inelastic scattering 
resulted from the inclusions of Au atoms in the indium- 
oxide matrix: By comparison, an undoped Y^O^-x sam- 
ple showed Tg" 1 that, at T w 0.5i"C, was more than 
three orders of magnitude smaller^. Since these 'pure' 
and the Au-doped samples had otherwise quite similar 
parameters (their and diffusion constant were the 
same to within 30% ), the non-trivial role of the gold 
in enhancing r^ 1 must be considered. It seems likely 
that this is a manifestation of the Sergeev-Mitin mech- 
anism for electron-phonon scattering in disordered met- 
als. The gold impurities in our samples are heavier than 
the host atoms, and being inert they are also loosely at- 
tached to the indium-oxide lattice. These factors limit 
their ability to follow the lattice movement, and thus the 
main assumption of the Sergeev-Mitin mechanism is ful- 
filled. At the same time, the Au atoms are active as 
local soft-modes which could be very effective in dephas- 
ing the electrons^. However, being weakly coupled to 
the lattice, they cannot efficiently dissipate the energy, 
gained by inelastic collisions with the electrons, to the 
bath. Therefore the Au inclusions contribute to dephas- 
ing much more than to energy relaxation. Note indeed 
that the phase coherence length in these samples is dom- 
inated by the interaction of the electrons with these local 
modest and it is ~ 0.4/im at T = 0.3K as compared with 
L ep w 100/xm. That the dephasing rate exceeds the en- 
ergy relaxation rate by many orders of magnitude is quite 
a general property of low temperature transport, which 
follows from the different temperature dependencies of 
energy relaxation processes on one hand and dephasing 
on the other hand. 

To summarize, we have employed a scaling analysis 
of non-ohmic resistance curves in order to extract the 
electron-phonon scattering rate of metallic films. The 
method makes use of the temperature dependence of the 
resistivity therefore it is best suited for those cases where 
the resistance can be used as a sensitive thermometer. 
Our analysis, which assumes quasi-two dimensional sam- 
ples may be easily extended to other dimensionalities. In 
these cases the temperature dependence of the resistance, 
at sufficiently low temperatures, is dominated by a power 
low behavior, R(T) = i?n(l — jT u ). From here it follows 
that the scaling function of long samples satisfies the 
relation: 



7 



1 



1 



Having the equilibrium parameters (7 and v), the exper- 
imental data of R(T, F) can be fitted to the above form 
and both rj and p can be extracted. The electron phonon 
length is then deduced from Eq. (|17|) . 

We reiterate that in order to apply the scaling ap- 
proach, the following conditions should be satisfied: (a) 
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FIG. 4: Left panel: The differential resistance versus the voltage V normalized by the temperature T using the data of Fig. 3. 
The full line is a best fit to the data using Eq. 1221 . Note the data collapse for the three lowest temperatures and the deviations 
from the short-sample behavior at T = 837m.K\ Right panel: Data from the same figure (Fig. 3) plotted as function of F(= V/L) 
normalized by T 2 . The full line is a best fit to the data using Eq. I20H . Note the data collapse for the two highest temperatures 
and the deviations from the long-sample behavior at T = 279mK. 



The heat transfer from the electrons to the bath is domi- 
nated by the electron-phonon collisions; (b) the electron- 
electron diffusion length should be much smaller than the 
energy relaxation length. 

On the other hand, the scaling approach is insensitive 
to the inclusion of other ingredients such as two level sys- 
tems and Kondo impurities as long as they do not serve 
as additional channels for heat conduction to the bath. 
Furthermore, for long samples the quality of the contacts 
is of minor importance, since the amount of heat trans- 
ferred by the electrons through the contacts is anyhow 



negligible. For short samples, however, our analysis as- 
sumes that the electrons near the contacts are at the bath 
temperature. This means that the contacts are ideal heat 
sinks, a caveat that should be borne in mind when using 
contacts made of a superconducting material. 
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